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Abstract 

I show how there is an ambiguity in how one treats auxiliary variables in gauge theories including general relativity 
cast as 3 + 1 geometrodynamics. Auxiliary variables may be treated pre-variationally as multiplier coordinates or as 
the velocities corresponding to cyclic coordinates. The latter treatment works through the physical meaninglessness of 
auxiliary variables' values applying also to the end points (or end spatial hypersurfaces) of the variation, so that these 
are free rather than fixed. [This is also known as variation with natural boundary conditions.] Further principles of 
dynamics workings such as Routhian reduction and the Dirac procedure are shown to have parallel counterparts for 
this new formalism. One advantage of the new scheme is that the corresponding actions are more manifestly relational. 
While the electric potential is usually regarded as a multiplier coordinate and Arnowitt, Deser and Misner have regarded 
the lapse and shift likewise, this paper's scheme considers new flux, instant and grid variables whose corresponding 
velocities are, respectively, the abovementioned previously used variables. This paper's way of thinking about gauge 
theory furthermore admits interesting generalizations, which shall be provided in a second paper. 
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1 Introduction 



Let (M, grA(x a )) be the spacetime topology and metric in the coordinates x a . Arnowitt, Deser and Misner (ADM) [T] 
([2] also has a good exposition) split spacetime with respect to a family of spatial hypersurfaces, T,\. Under this split, x a 
splits into (A, x^), where A a coordinate label time for each spatial hypersurface in the split and x^(X) a coordinatization 
of that spatial hypersurface. Then ADM take grA to split as follows: 

grA = /? 7 h lS J ■ (1) 

Here, h~fg(\, x^) is the induced metric, i.e. the metric induced by the spacetime metric grA on the spacelike A-hypersurface. 
a(A, x^) is the lapse, which relates the passage of coordinate time between infinitesimally adjacent surfaces T,\ and S^+dA 
to the passage of proper time r: 

dr = a{\,x~*)d\ . (2) 
/3(A, x M ) is the shift, which relates the coordinates on E>,+dA to those on Ea: 

x 7 (A + dA) -x 7 (A) = -/? 7 dA . (3) 

I.e., the lapse is the change in spacing in how the spatial hypersurfaces are stacked, while the shift is the displacement 
made in identifying the spatial coordinates of one slice with those on an adjacent slice. Thus these two can be envisaged 
together as the determiners of what extrinsic geometry each spatial hypersurface T,\ is to have within spacetime (see e.g. 
p 346 of [3])- The ADM split is then applied to the spacetime form of the general relativity (GR) action to produce an 
ADM action for split space-time|]] 

IadmI^,^,/^,^] = y*dA J d 3 xL(/i / _ 11/ ,a;/i AI „,/3 M ] = ^dA J d 3 xVha j — ^— ^ £pK„}{h pcr £ph pa } + ^ 

(4) 

The ADM split and variants have important applications in quantum gravity (see e.g. [4j [6l [3 [8] ) and in the study of the 
compact binary problem in numerical general relativity (see e.g. [5]). 

It is 'standard lore' in varying the ADM action to regard the lapse and the shift as coordinates; as the corresponding 
velocities do not feature in the action, these are Lagrange multiplier coordinates. Denoting Lagrange multiplier coordinates 
in general by m& , variation with respect to these produces simplified Euler-Lagrange equations of the form 

8L 0. (5, 



dniB 



This is how the Hamiltonian and momentum constraints of GR arise in the ADM formulation (see the Appendix) . 

Now, it is a fairly standard observation (see e.g. [10]) that the above expressions for the lapse and shift can, at the 
post-variational level, be straightforwardly rearranged to be interpreted as velocities, 

a(\,x-y) = ^- , /T(A,^) = -^ . (6) 

What's new in this paper is to regard the objects of which the lapse and the shift are the velocities as new variables 
prior to variation. Reading off (|44[) . these are the proper time r attributed to the spatial slice in question and the grid 
of coordinates x^ on that slice, but in their capacity as new variables to be varied I term these the instant /(A,:^) and 
the grid F r (X,x^ 1 ). [I use an 'i 7 " to denote this as it is a subcase of a more widely applicable concept of spatial fi 



rame 



variable]. Through showing this to work out acceptably as an alternative interpretation, this paper shows that there is 
in fact an ambiguity in the long-standing, well-known and much-used statement that the lapse and shift are Lagrange 
multipliers 

Moreover, the alternative better complies with relational first principles, which is of interest along the lines of estab- 
lishing GR to have multiple foundations, along the lines of Wheeler's arguments for 'many routes to GR' 01 mill The 
relational perspective, implements ideas of Leibniz [12] and Mach [13j to modern physics along the lines of Barbour's work 
[TJJ [16] ■ One starts with a configuration space Q of (models of) whole-universe systems. 



1 I use ( ) for function dependence, [ ] for functional dependence, ( ; ] for a mixture of function dependence before the semi-colon and 
functional dependence after it. I use lower-case Greek letters for spatial indices, except that a and /3 are reserved to denote lapse and shift. 
h = h(h ll ^(x v ')), D M , and R(x n ; h^(x w )] are, respectively, the determinant, covariant derivative and Ricci scalar associated with the 3-mctric 
h t j, v (x 7T ). The dot denotes d/dX from Sec 4 onwards. £p is the Lie derivative with respect to /3 7 . G a /3jS = rr tiv^ pa o ^fJ-vfopd } ^ th.6 

DeWitt supermetric [4] with inverse 

G ^pa _ ^Th^pfrva _ f l ^h'" T } which plays the role of kinetic metric (see e.g. \5\) in GR. Overbar 
denotes densitization (multiplication by \fh) and underbar denotes dedensitization (division by \fh). 

2 For, different routes offer different techniques and different insights, suggest different alternatives or generalizations against which the 
established theory can be tested, and different routes to a given classical theory may be quantum-mechanically inequivalent. 
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Then configurational relationalism is the declaration that certain motions acting on Q are to be physically meaningless. 
One wajffl of implementing this is to use arbitrary-G-frame-corrected quantities rather than bare Q configurations, where 
G is the group of physically meaningless motions. Despite this augmenting Q to Q x G, variation with respect to each 
adjoined independent auxiliary G-variable produces a constraint which wipes out one G variable and one redundancy 
among the Q variables, so that one ends up on the quotient space Q/G (the desired reduced configuration space). This 
is widely a necessity in theoretical physics through working on the various reduced spaces directly often being technically 
unmanageable. 

Temporal relationalism is the notion that there is no meaningful primary notion of time for the universe as a whole. 
One implementation of temporal relationalism in modern physics is through using manifestly reparametrization invariant 
actions, though Barbour would prefer this to be attained without any extraneous time variables either. If one is working 
towards deriving GR in this picture, one would e.g. take Q to be the space of positive-definite 3-mctrics on a fixed 
spatial topology, Riem(E), and G to be the 3-diffeomorphisms Diff(E) that correspond to the coordinatization of space 
being physically meaningless. It is then clear upon inspecting the split space-time GR action that the instant and grid 
interpretations of the auxiliaries therein give the action a manifest reparametrization invariance while the lapse and shift 
multiplier coordinate interpretations obstruct this by inhomogenizing the Lagrangian's dependence on d/d\. To have 
no extraneous time variables either, one would require not an ADM-type action with its manifest lapse or velocity of 
the instant, but a Baierlein-Sharp-Wheeler (BSW) [T!5] or BFO-A0 type action in which this has been eliminated (or, 
alternatively, such an action could be adopted as one's starting-point). This paper completes this work by mapping out 
how the instant can be eliminated from the instant-grid version of the split GR action so as to produce the BFO-A action. 

There is an obstruction that explains why this kind of alternative interpretation of auxiliary variables has not been 
spotted until recently. At first sight, it looks preposterous to treat the lapse and shift as velocities prior to variation. This 
is because cyclic coodinates cb are well-known to lead to simplified Euler-Lagrange equations of the forn0 

- — = G B , constant , (7) 
oc B 

which are generally not equivalent to the ([5]) produced by the preceding formalism. However, there is a subtlety. In 
a gauge formulation of a physical theory, the configuration variables contain both dynamical degrees of freedom and 
nondynamical/auxiliary variables. Then by the gauge principle (or configurational relationalism), the reference frame 
variables are to be taken to be unphysical, i.e. arbitrarily specifiable at each instant, with no measurable quantities 
depending on which choice of these is made. This includes among the things that are to be arbitrary the values of the 
auxiliary arbitrary G-frame variables at the end points (or end spatial hypersurfaces in the case of field theory) of the 
variation. Thus it is free end point variation (also known as variation with natural boundary conditions) |22| that turns 
out to be the correct type of variation as regards reference frame variables in gauge theories. In Sec 2, I show that this in 
no way affects the outcome of varying with respect to multiplier variables, but it does alter the outcome of varying with 
respect to cyclic variables precisely so as to remove the above inequivalence. 

In Sec 3, I then show that there is also no difference between multiplier elimination in the auxiliary multiplier formula- 
tion and cyclic velocity elimination (by the more subtle procedure of Routhian reduction, see e.g. [5]) in the above auxilary 
cyclic coordinate formulation. I also explain how one subsequently gets an alternative to Dirac's appending procedure 
[23] . These sections are supported by relational particle mechanics toy models [TH [15l [24l [25] [26] [27] - zero total angular 
momentum universes formulated by frame-correcting with respect to the rotations with a rotational auxiliary analogue 
of the grid-shift object. Armed with all this, in Sec 4 I give the manifestly reparametrization-invariant alternative to the 
ADM formulation of Einstein-Maxwell theory. This requires, in addition to the lapse and shift being re-interpreted pre- 
variationally as instant and grid variables, for the conventional electric potential $ to be re-interpretated as the velocity 
of a new (corrected magnetic) flux variable ^ (which is another example of frame variable). Then performing Routhian 
reduction on this action so as to eliminate the instant leads to the recovery of the BFO-A action that complies with the 
first principles of the relational approach, which working is in direct parallel to how multiplier elimination of the lapse 
from the ADM action gives the BSW action. I conclude in Section 5 with a discussion of applications and extensions, 
among which the case of auxiliaries more general than multipliers or cyclic coordinates is presented in a subsequent paper. 
The Appendix provides Einstein-Maxwell theory's ADM and BSW formulations (and the BSW elimination that links 
them) for useful comparison. 

2 An alternative picture for gauge formulations 

One common way of formulating gauge theories is in the Lagrangian picture, in which one builds Lagrangians that possess 
certain groups of symmetries, G. Another common way of formulating gauge theories is in the Hamiltonian picture, in 
which one appends certain constraints (e.g. in the sense of Dirac [53]) that are related to the action of the gauge group. 
The Dirac procedure has the additional advantage that even if one does not know what symmetries the theory is to 

3 This is along the lines of |17l . which is closely connected with well-known approaches to gauge theory. For Barbour's own and somewhat 
conceptually different way of thinking about configurational relationalism ('best matching'), see e.g. [14 18 . 
4 This action l|44j l first appeared in |17| . building on |16| and insights in [20] , See also [21) for discussion. 
5 In Sees 1-3 the dot denotes d/dA. 
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possess, the procedure nevertheless systematically derives the theory's constraints (see e.g. [2E])- One can move between 
the two pictures by Legendre transformation. One approach within the Lagrangian picture is to attempt to impose a 
group G of symmetries on a Lagrangian, and then use the Dirac procedure to systematically determine whether this 
imposition is consistent and completed 

One way of building a formulation in which the group of symmetries G (whether actual, expected or imposed) is 
present and manifest is to enlarge the configuration space Q to Q x G by adjoining auxiliary variables g&- These show 
up in the Lagrangian formalism as arbitrary G-frame corrections to the Q- variables [17[ 121] . To see how this gives gauge 
theories, consider for the moment that the auxiliary variables take their usual guise as multipliers in the generalized sense 
(coordinates whose velocities do not occur in the action). E.g. one conventionally views thus the electric potential 
and the shift /3 7 and lapse a of GR split with respect to spatial hypersurfaces. Further examples are the translation 
and rotation auxiliaries, a and b in mechanics models that are spatially relational with respect to the Euclidean group 
[HI [25l [26] (as opposed to being based on the supposition of absolute space). 

Variation with repect to the auxiliaries <?b then produces 

which are the constraints that one is expecting if one is indeed to obtain thus a gauge theory with gauge group G. These 
constraints then use up both the auxiliary variables' degrees of freedom and an equal amount of degrees of freedom from 
the Q, thus indeed leaving one with a theory in which the dynamical variables pertain to the quotient space Q/G (i.e. 
variables among the original Q which are entirely unaffected by which choice of G-frame is made). 

As a simple example, consider the translation and rotation invariant mechanics p31 [25l [26] that follows from the 
Jacobi-type action 

\[q v q v a,b] = J d\L(q I ,q I ,a,b)=2 J d\ ^T{V + E} 



with kinetic term 



potential term 



(9) 

T{gj,gj,a,b) = -^milq! - a-b x qj 2 , (10) 



i=i 



V = V(|r/. 7 |alone) (11) 

and total energy E. Then variation with respect to a gives the multiplier equation = dL/da = Yli^iPji i- e - zero total 

momentum, and variation with respect to b gives the multiplier equation = dL/db = ^¥=i1t x Pri l - e - zer0 total 
angular momentum. See the Appendix for another example (the standard ADM split of Einstein-Maxwell theory) . 

The main point of the present Paper, however, is that taking a, b, <£>, /3 7 and a to be multipliers is not the only viable 
interpretation. They may also be interpreted as velocities associated with an auxiliary cyclic coordinate. This is possible 
because the variables in question are auxiliary. Thus by the gauge principle, their values are entirely arbitrary. Thus, in 
particular, they take arbitrary rather than fixed values at the end points of the variation. Thus the appropriate type of 
variation is free end point (FEP) variation (also known as variation with natural boundary conditions) [22j . FEP variation 
involves more freedom than standard variation, i.e. it involves a larger space of varied curves, 

q A (X,/j,,u) = g A (A) +fjY(X) + vZ(X) such that Y(X l ) = G and Z(X f ) = (12) 

(where fj,, v are further parameters). Then consider configurations {<7a,3b} for gB entirely arbitrary through being 
unphysical. With this in mind, I start from scratch, with a more general Lagrangian in which both the auxiliaries gB and 
their velocities gs in general occur: L(A, 9a, 5b, 9a, <7b)- 
Then variation with respect to <?b gives 

f Xf f Xf ( dL d f dL 1 1 f dL 1 A/ 

= 51=/ dXSL(X,qA,g B ,q A ,g B )= / dX6g B \ ^- - — ^ + ^ffB ■ (13) 



\i J\> I c*9b dX I dgB J J I og B J A 

Because gB is auxiliary, this variation is free, so Sgs |a ; and SgB\\ f are not controllable. Thus one obtains 3 conditions per 
variation, 



or, in terms of momenta p E 



dL 


d 


(dL] 


dL 


dgB 


~ dX 




dgB 


dL/dg B , 










dL 








dgB 


= P B , 


P B k = 



dgB 



(14) 



p B L ■ (15) 



6 Constraints that do not correspond to G could arise as integrabilities for other constraints in the theory, signalling that G need be extended 
if one is to get a consistent theory. See the end of Sec 4 for an example of this. So many constraints can arise in this way that the resulting 
"theory" ends up trivial (all degrees of freedom used up) or inconsistent (more independent constraints than degrees of freedom) |16| . 
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If the auxiliaries g& are multipliers me, (|13p reduces to 



dl , . 

PB = , = 16 

and redundant equations. I.e. the A^ and Xf terms automatically vanish in this case by applying the multiplier equation 
to the first factor of each. This is the case regardless of whether the multiplier is not auxiliary and thus standardly varied, 
or auxiliary and thus FEP varied, as this difference translates to whether or not the cofactors of the above zero factors 
are themselves zero or not. Thus the FEP subtlety in no way affects the outcome in the multiplier coordinate case. 
If the auxiliaries gs are cyclic coordinates cb, the above reduces to 

P\=° = P% (!7) 

and 

p B = , (18) 

which implies that 

P B = C , (19) 

but C is identified as at either of the two end points (fTT|) . and, being constant, is therefore zero everywhere. Thus (|18p 
and the definition of momentum give 

J^ B = °- < 2 °> 

So, after all, one does get equations that are equivalent to the multiplier equations ([5]) at the classical level, albeit there 
are conceptual and foundational reasons to favour the latter, as laid out in the Introduction. For previous discussion of 
this case in the literature, see [2T1 |2"01 [T71I5T] . 

For the translation and rotation invariant mechanics example, this works out as follows. The appropriate action is 

\[q v q_ v aM = J d\L(q I ,q J ,a,b) = 2 J d\^/T{U + E} , (21) 

with 

1 N 

~t{g v QjAS = ^^2 m i\<ii -a-b x £ 7 | 2 (22) 

1 7=1 

and U = —V. Then variation with respect to a gives C = dl/da = Ei=i2j and the FEP conditions dL/da\ Xt = = 
dL/da\x f , so C — dL/da = at the end point, but C is constant, so C = everywhere, so one obtains Y^,i=iPj = 
again. And variation with respect to b gives D = dL/db = J2i=i Qj x Pj an d FEP conditions dL/db\\ i = = dL/db\\ f , 
so D = dL/db = at the end point, but D is constant, so D — everywhere, so one obtains E7L1 Qj x P 7 = again. 
For further examples of this kind of approach to gauge theories, see e.g. Sec. 4 (3 + 1 split GR), Sec. 5 (3 + 1 split of 
Einstein-Maxwell theory), [271 126] (mechanics that is scale-invariant as well as translation and rotation invariant), and 
[32] [20], 131] alongside Paper II for a conformal theory of gravity and conformogeometrodynamical formulations of GR 
(conformal split of 3 + 1 split GR along the lines of York's work [33] on the initial- value formulation). 

The case of ([14"]) in which the auxiliary is neither cyclic nor a multiplier I leave for Paper II [31], as the confor- 
mogeometrodynamical formulations in which it applies [20, 31] are rather more complicated than the present paper's 
examples. 

As regards how one determines in the first place which variables are the auxiliaries that one is to FEP vary with 
respect to, it turns out that for the type of examples in this suffices to FEP vary with respect to the gs, as follows. While 
some sets of variables contain partly physically relevant and partly physically irrelevant degrees of freedom, the gauge 
principle only gives license to vary in the FEP way with respect to sets of variables that are purely physically irrelevant. 
So FEP variation should be applied only to isolated physically irrelevant variables. Moreover, one does not have to isolate 
all of these, as follows. In gauge theory, physically irrelevant variables come in pairs that are associated with a single 
constraint, which, when taken onto account, ensures that neither member of the pair remains in the theory's equations. 
Variation with respect to either member of the pair produces the same constraint. Thus to get all of the constraints, one 
needs to isolate a full half-set of auxiliary variables i.e. to isolate one representative from each of the above pairs. But 
in the arbitrary G-frame method as used in this paper, the frame variables that one adjoins are clearly both such a 
half-set and already-separate. Thus for the type of theories considered in this paper, it suffices to FEP vary with respect 
to the gs- 



3 Establishing further equivalences between auxiliary cyclic coordinates 
and auxiliary multiplier coordinates 

We also need to establish that the a priori distinct procedures of cyclic velocity elimination (known as Routhian reduction) 
and muliplier elimination are equivalent. Consider L(qA, <Za, <?b) where gs are auxiliary coordinates. If these are taken 
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to be multiplier coordinates tob, then variation yields = SL/Stob. If this is soluble for tob, one can replace it by 
nig = tob(<Za,Qa), and then substitute that into L. If the gs are taken to be velocities corresponding to cyclic coordinates 
cb, then FEP variation yields 

= |^ = p B = C B . (23) 

OCb 

This is soluble for cb iff the above is soluble for tub- However one now requires passage to the Routhian in eliminating cb 
from the action: R(<7a,<7a) = L(gA, <Za, ca) — cbp B - But the last term is zero, by (|2"3"j). so Routhian reduction in the cyclic 
velocity interpretation is equivalent to multiplier elimination in the multiplier coordinate interpretation. 

Next, I consider how the Dirac procedure implies encodement by auxiliaries at the level of the Lagrangian. From the 
bare (G-uncorrected Lagrangian) L bare (<?A, 9a), Legendre transformation gives H(<7a,.P A ) = 4aP A — L barc (gA, q A ) for which 
Cb = are discovered by the first part of the Dirac procedure. Then, in the multiplier interpretation, 

H(q A ,p A , tub) = q*p A - L barc (g A , q A ) + m B C B (24) 

so that this procedure standardly follows from the standard Dirac procedure of appending constraints multiplied by 
Lagrange multipliers. In the cyclic velocity interpretation however, one has a new alternative to Dirac appending that 
involves appending constraints multiplied by velocities of cyclic coordinates. Then one is not forming a Hamiltonian, but 
rather an 'almost-Hamiltonian', in that it depends not only on positions and momenta but also on the velocities of the 
auxiliary cyclic coordinates, 

A(g A ,p A ,c B ) = q AP A + c B C B - L w (g A ,g A ) • (25) 

There is no problem uplifting Sec 2 or 3 to field theory. [Now one has free end spatial hypersurfaces (FESH) and 
C(x IM ) in place of C, constant, but all the above arguments straightforwardly carry through]. 



4 Manifestly reparametrization invariant alternatives to the ADM and 
BSW-type formulations of Einstein— Maxwell theory 

As delineated in the Introduction and Appendix, in the ADM split, lapse and shift variables are used and have the status 
of Lagrange multipliers. In this paper, instead, I use the new split 

'P p pp - i 2 p v 

//!. a \ K C I *•.><; i 

where is the frame (grid) and / is the instant, the meanings of which are explained in the Introduction. The consid- 
erations in Sec 2 shall ensure that the subsequent variation works out. 
The GR kinetic term is now 

T = G^ pa { V - £ph ia ,}{h pa - £ph pa } . (27) 

and the GR potential term is R as before. 

For greater generality, I include the electromagnetic field (massless spin 1 field) as an example of matter fieldQ In 
Maxwell electrodynamics, the usually-employed electric potential variable is defined by 

V$ = -E-A. (28) 



Here A is the magnetic potential such that B_ = V x A, where B_ is the magnetic field and E_ is the electric field. The 

r(t) 



variable I use in this paper is, rather, $ such that 4/ — <I>. In the same way that $ is a 'work per unit charge' J 3 



\& is an 'action per unit charge' J dt J r ^ — ■ Another interpretation is as follows. In the flat spacetime context, 
integrating (|28[) and applying the Faraday-Lenz law to the first factor and Stokes's theorem, mixed partial equality and 
the definition of A to the second factor, 



a r r r r /a corrected rate of change N * 

— / / B-dS- I 
r(t) dt J J s{t) J J S (q 



§ = -j) {Eft + A^jdx^ = — I I B-dS- I I B-dS=\ of magnetic flux ) (29) 



(for surface S with boundary curve T). By Leibniz's rule and theorem this is also 



/ change in magnetic flux due to 
$ = / / B ■ dS_ = change in shape of the surface over time ] . (30) 



as(t) 



7 Including scalar fields is straightforward, and inclusion of spin-1/2 fermions follows along standard lines. This paper's treatment of 
electomagnetism generalizes to Yang-Mills theory and to the various associated scalar and spin-1/2 gauge theories. 
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Thus 



*=// B-dS- dt / B ■ dS = (corrected magnetic flux) , (31) 
J Js(t) J J Js(t) 



IS(t) J J JS(t) 

the first term being magnetic flux and the second term being the correction; this is also 

magnetic flux due to change in shape 



/r r i magnetic nux auc to cnangc m snapc \ 

dt / / ^ B_ ■ dS_ = I of the surface over time integrated over time J . (32) 

Electromagnetism in this formulation has kinetic and potential terms 

T(A,M = \\A,-dM 2 ■ y[A^] = \F^ (33) 

for F^y = d^Ay — dyA^ the spatial part of the electromagnetic field strength tensor. This formulation does indeed work in 
the flat spacetime context: C(x fl ) = SL/S'i = c^IT 1 but C on the end spatial hypersurface is and C is time independent, 
so C is everywhere, and thus we recover the Gauss constraint, d^H^ = 0. 

I then obtain the suitable manifestly rcparamctrization invariant action for Einstein-Maxwell theory by building the 
above up to hold on curved spaces and then coupling it to GR so that it holds on dynamically changing curved spaces. 
The complete Einstein-Maxwell action is then 

l A [V,V,^,i,-F„i,f] = J d\J d 3 ar/L A (V,iV; K^A^F^] = j d\ j d 3 xVhi + ^ + R - V A J 

fdx fd 3 ^i f z ^VHV - £ F h P «} h ^iK QjMiK -£pA v - d^} 

J J X \ AP 2P 

-\f^F^ + R^ (34) 
From this starting point, the conjugate momenta are 

^ = SU = yft Glivpa{h _ £ h } and n ^ dl = Vh {A ^ _ £ ^ _ (35) 
dh^ 21 dA M I 

Variation with respect to / gives as a secondary constraint 

C{x?) = % = G^ir^ir^ + -^f I^IF + ^F^F^ - VhR , (36) 
ol 2\Jh 4 

and the FESH conditions dLp^/dl^ = = d\-A/dI\\ f , so C(x^) = dL^/dl = at the end spatial hypersurface, but C is 
hypersurface-independent, so C = everywhere. Thus one arrives at 

G^tt^tt^ + ^n M IF + ^F^F^ -VhR = Q, (37) 

i.e. the Einstein-Maxwell Hamiltonian constraint, H A . Variation with respect to ^ gives 

G{x*) = ^ = D^W = d^W , (38) 
aw 

and the FESH conditions SLa/c^Ia, = = <9La/<9 , J , |a / , so G(x fl ) — dLx/dif — at the end spatial hypersurface, but G 
is hypersurface-independent, so G = everywhere. Thus one arrives at 

D^W = , (39) 

which is the electromagnetic Gauss constraint, Q, of Einstein-Maxwell theory. Variation with respect to F^ gives (modulo 
a Gauss constraint term) also as a secondary constraint 

E 7 {xn = |^ = -2£„V - n s {D s A, ( - D^A S } (40) 

and the FESH conditions <9La/<9F 7 |a, = = d\-x/dFi\\ f , so E J (x fJ -) — <9La/<9F 7 = at the end spatial hypersurface, but 
E 1 is hypersurface-independent, so E 1 = everywhere. Thus one arrives at the momentum constraint of Einstein-Maxwell 
theory, U^, 

- 2D s ir 7 5 = li s {D s A 1 - DjAs} . (41) 
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This formulation's equations of motion t\^ v = <9La/ dh pu , n p = dLj^/dA^ propagate the above constraints without giving 
rise to any further ones. 

Taking the Lagrangian form of (|36|) , one can now solve for /: 



i = ±lJl+lA (42 ) 
2V R+Ua 

(and take the + sign by convention of direction of forward march of time) , from which / is straightforwardly algebraically 
eliminable from the action (fS"6|) by Routhian reduction. Thus one recovers the 'BFO-A' action 



Ibfo-a[V/> h^A^A^F^] = J dxj d 3 sL BF o-A(V,i M ; h^^^F^} = J dAd 3 a;- v /{T + T A }{R + U A } = 

(43) 



dA J d 3 xVh 



\ 



2/i^{ i„ -^-^{i-^-^f} i?- I . (44) 



This is my formulation's counterpart and equivalent of BSW's multiplier elimination procedure [19j . thus revealing [35] 's 
"explicit use is made of multiplier status " to be an alternative rather than an obligation. 
Taking this action as starting point, the momenta are 



- = = Ji±^G^{ V - £ F h, v } , W = = VhJi±^h^{A v - £ F A V - dM (45) 



9V V T + Ta " ^ J '1T + T, 

now one does not have not a priori an instant / of which the rate appears in the equations, but rather an emergent 

1 / R+V 



quantity ^ * / j+Ta ' Then there is a primary constraint 



+2hrWW =h^hMypj^{A p - £ p A p ~ d p i>}h^^±^{A a -£pA a - dA}= V^±^{T+T A } = Vh{R+U A } 

. (46) 

(by the definition of momentum, and using that G Q( a 7 5 is the inverse of G a ^ s ), the Hamiltonian constraint. Variation 
with respect to *]/ gives 

G(x») = ° B f T °- A = D p W = d^W , (47) 
aw 

and the FESH conditions SLbfo-a/ = = SLbfo-a/ 'd^\\ f , so G(x M ) = 9Lbfo-a/5 v E' = at the end spatial 
hypersurface, but G is hypersurface-independent, so G = everywhere. Variation with respect to F p gives (modulo a 
Gauss constraint term) the secondary constraint 

E 7 {x») = B J°~ A = -2iW« - n A '{D 5 A 7 - D 7 A 5 } (48) 

and the FESH conditions <9Lbfo-a/<9-F' 7 |a 1 = = d\-BFO-A/dF~«\ Xf , so E 1 {x^) — <9L B fo-a/<9F 7 = at the end spatial 
hypersurface, but Fj is hypersurface-independent, so F^ — everywhere. Thus one arrives at 

- 2A,tt/ = n A '{^A 7 - D 7 A 5 } , (49) 

the momentum constraint of GR, 7i M . The BFO-A equations of motion ■k^ 11 ' = (9Lbfo-a/c?/i m „ propagate these constraints 
without giving rise to any further ones. 

Thus in this scheme temporal relationalism as implemented by manifest reparametrization invariance without ex- 
traneous variables gives the Hamiltonian constraint of GR and spatial relationalism as regards 3-diffeomorphisms as 
implemented by arbitrary G-frame cyclic coordinate velocity corrections to the metric velocities gives the GR momentum 
constraint. In fact, even if a number of features of Lbfo-a are not known, if it is built according to temporal rela- 
tionalism, the Dirac procedure is applied and consistency is demanded, Ibfo-a emerges as one of very few possibilities 
[TBI |3"S1 [T71 [2Tj . The above very clear way of obtaining the GR constraints (and lack of prior knowledge of the form of 
the GR equations to do so) and the accommodability into this scheme of a sufficiently wide range of fundamental matter 
fields to describe nature at the classical level [Ml IIZI HI] make 1(13]) itself interesting as a starting point for gravitational 
physics. In fact, starting with temporal relationalism (and a configuration space of spatial 3-metrics) alone suffices [50] , 
as that produces 7i as above, which then happens to dictate TC and Q as integrability conditions (see also [29]), whose 
subsequent appending or encoding leads to the recovery of ([44]) Fl Though this may be considered to involve luck, as, in 
other cases [3D], relying on integrability can cause constraints to get 'missed out'. 



"Strictly speaking, these references consider the BSW formulation counterpart of these results; the BFO-A counterpart to these results is a 
new result of the present paper. 
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5 Conclusion 



In this paper, I have explained that variation with respect to auxiliary variables is free end point (FEP) in finite cases 
or free end spatial hypersurface (FESH) in field-theoretic cases. This makes a difference to the outcome of varying cyclic 
coordinates if they are auxiliary, whereby interpreting an auxiliary present in one's action as a cyclic velocity is equivalent 
to interpreting it as a multiplier coordinate Thus one has an alternative to the ADM split: instead of lapse a and shift 
0ft multiplier coordinates, one has instant / and grid cyclic coordinates, which are related to the a and /3 M by I = a 
and Fft = fl^. Also one has an alternative to the usual electric-magnetic split formulation: instead of an electric potential 
multiplier coordinate one has a flux cyclic coordinate \& that is related to it by \E r = <&. 

Further principles of dynamics manipulations continue to give the same answers in the new scheme. For example, 
Routhian reduction now works out to be equivalent to what was previously multiplier elimination. E.g. the new scheme 
has a Routhian reduction instant- velocity elimination counterpart of Baierlein, Sharp and Wheeler's multiplier elimination 
of the lapse. As another example, this scheme supplants Hamiltonians H by 'almost Hamiltonians' A which depend on 
positions and momenta for the usual variables but on postions and velocities for manifestly pure- frame variables. There 
is then no difference in constraints arising from the multiplier picture's H and the cyclic coordinate picture's A in the 
position-momentum language for the theories in this paper (for which the auxiliary variables, however represented, do 
not enter the position-momentum form of the constraints). Nor is there any difference in the thin sandwich conjecture 
(solving the Lagrangian form of H.^ with matter struck out for /3 M |38j or now for i 7 ^, or the Lagrangian form of TL^ and 
Q for (3^ and $ [35] or now for and 

One advantage of this paper's formulation is that it is more manifestly temporally relational than ADM's. It is true 
that the instant-less BFO-A formalism is 'even more' manifestly temporally relational in Barbour's sense of not containing 
an extraneous time variable, but, on the ther hand, thinking along the lines of this paper's formalism and then eliminating 
the extraneous time variable permits the relational program to go beyond its current spatially compact without boundary 
setting to cases with boundaries or that are open. 

The cyclic velocity picture gives a different account of how the frozen formalism problem of GR (see e.g. [71 IBl 1571 124) ) 
arises. 

w = S = -i s -n' (50) 

dl dl 

by, respectively, trivial direct computation, the chain rule applied to the Legendre transformation relation and the defi- 
nition of momentum. Thus 

Tl I +H = Q. (51) 



Similarly, 

Il£ + Hft = . (52) 

At this stage, these contain linear momenta, so, in particular, the first quadratic equation has no trace of the frozen 
formalism. It would amount to 

ih^f=H\^> (53) 

at the quantum level, which is an (instant, or proper,) time-dependent Schrodinger equation rather than a stationary 
equation like the Whcclcr-DcWitt equation. However, one has also the FESH condition n 1 = 0, which precisely wipes 
out the linear momentum in the crucial, quadratic constraint, leaving 

H = 0. (54) 

At the quantum level this amounts to the usual stationary Wheeler-DeWitt equation 

= . (55) 

However, in cases/simplified models in which there is a privileged background, the situation here is, very transparently, 
that FEP is then inappropriate, whereby some kind of time-dependent Schrodinger equation survives. As many quantum 
schemes proceed via the position-momentum form of the constraints that is unchanged by passing to this paper's for- 
malism, these schemes are likewise unchanged. I leave investigation of whether this article's change has any effect on any 
other quantization procedures for a future occasion. 

Finally, I explained how auxiliary variables can be more general than either multiplier coordinates or the velocities 
corresponding to cyclic coordinates, which affords further extension of how gauge theories are treated in this paper. This 
taking one further afield than standard formulations of standard theories of physics (e.g. these more general cases apply 
to [201 131]), I present these in a second paper [53] . 
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Appendix: ADM and BSW-type formalisms of Einstein— Maxwell theory 

The conventional 3+1 split action for Einstein-Maxwell theory is 

Iadm[V> h^^^A^P^a,®] = J d\ J d 3 a;LA(V, a \ K V ,A^ /3 M ,$\ = J d\ J d 3 xVha | - + J A +R - ^F^F^A 

J J X a \ 4a 2 2a 2 

-■^ + ^} ■ (56) 

Here, a is the lapse, /3 M is the shift and <!> is the electric potential; these are all taken here to be Lagrange multipliers. 
See footnote 7 for the rest of the notation. 
From this, the conjugate momenta are 

tt^ = = * G »»P°{} 1 - £ h } an d IP = A- = — /^{^ - jfi^ - d,<P} . (57) 

2a OA,, a 

Variation with respect to a gives the Einstein-Maxwell Hamiltonian constraint, Tt A : 

G^-k^-k"" + ^n^IP + ^F^F^ - VhR = . (58) 

Variation with respect to $ gives the Gauss constraint, Q: 

Z^IP = . (59) 

Variation with respect to F^ gives (modulo a Gauss constraint term) the Einstein-Maxwell momentum constraint Tt A : 

- 2D S 7T 7 S = n s {D 5 A 7 - D 1 A S } . (60) 

One then finds that propagation of these by the ADM equations of motion for Einstein-Maxwell theory, tt^" = 9Lbsw /dh^, 
IE' = OLbsw/QA^, gives no further constraints. 

Using the Lagrangian variables version of the a-multiplier equation, 



(and taking the + sign by convention of direction of forward march of time), a is straightforwardly algebraically eliminable 
from the ADM action giving the Baierlein-Sharp- Wheeler (BSW)-type action for Einstein-Maxwell theory: 

Ap, , A^ , , j d\ J d 3 xLBsw (V, A„; h^,A^,^, $] = J dAd 3 .xV{T + T A }{i?+U A } = (62) 



dA J d 3 xVh 

From this starting point, the conjugate momenta are then 



\ 



2h^{A ll - £ p A li - d^}{A v - £ P A U - )\R- -F^F^ \ . (63) 



= ^ = >/^G'-{^ - £M and IF ^ = VhJ^±^h^{A v - £ A V - £„$} . (64) 
dh^v V I + I a oA^ V I + I a 

Then there arises as a primary constraint the Hamiltonian constraint of Einstein-Maxwell theory (|58p . Variation with 
respect to $ gives the Gauss constraint (|59"1) . Variation with respect to F^ gives the Einstein-Maxwell momentum 
constraint (|60[) (again modulo a Gauss constraint term). One then finds that propagation of these by the BSW equations 
of motion 7T MI/ = <9Lbsw/<9/i m „, IP 1 = 9Lbsw/9A m , again gives no further constraints. 
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